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We study theoretical bases of the dual Higgs theory for confinement physics in QCD 
in terms of monopoles and the gluon configuration in the maximally abclian (MA) 
gauge. Abelian dominance for the confinement force can be analytically proved 
by regarding the off-diagonal angle variable as a random variable in the lattice 
formalism. In the long-distance scale, the contribution of off-diagonal gluons to 
the Wilson loop cancels each other and exhibits a perimeter law behavior, which 
leads to exact abelian dominance on the string tension if the finite size effect of 
the Wilson loop is removed. We investigate the appearance of the monopole in the 
QCD vacuum, considering the role of off-diagonal gluons. The monopole carries a 
large fluctuation of the gluon field and provides a large abelian action in abelian 
projected QCD. Due to the partial cancellation between the abelian part and the 
off-diagonal part of the QCD action, the monopole can appear in QCD without 
large cost of the QCD action. The off-diagonal gluon is necessary for existence 
of the monopole at the short-distance scale. We study monopole condensation, 
which is the requirement of the dual Higgs theory, by comparing the QCD vacuum 
with the monopole-current system. We find that "entropy" of monopole-current 
dominates than its "energy", and the monopole seems to be condensed at the 
infrared scale in the QCD vacuum. 

1 Dual Superconductor Picture for Confinement in QCD 

The strong interaction is subjected to Quantum Chromodynamics (QCD). Due 
to the self-interaction of the gauge field, the QCD gauge coupling becomes very 
strong in the low-energy region, while it is weak in the high-energy region. 
Accordingly, QCD phenomena are divided into two theoretical categories: the 
strong coupling leads to complicated nonperturbativc phenomena such as color 
confinement and chiral symmetry breaking, while high-energy phenomena are 
understood by the perturbative QCD. Since there is no well-established ana- 
lytical method for nonperturbative phenomena, one must carry out the Monte 
Carlo simulation based on the lattice QCD or apply the effective models de- 
scribed by the relevant degrees of freedom for the low-energy physics. As for 
the chiral dynamics, the pion and the sigma meson, which are composite modes 
of quark and anti-quark, play an important role for the infrared effective theory 



such as the (non-)linear sigma model and the Nambu-Jona-Lasinio model. On 
the other hand, confinement is essentially described by the dynamics of gluons 
rather than quarks. 

In 1970's, Nambu proposed an interesting idea that quark xonfincmcnt 
can be interpreted using the dual version of the superconductivit5a. In the su- 
perconductor, Cooper-pair condensation leads to the Meissner effect, and the 
magnetic flux is squeezed like a quasi-one-dimensional tube as the Abrikosov 
vortex. On the other hand, from the Regge trajectory of hadrons and the lat- 
tice QCD, the confinement force between the color-electric charge is brought 
by one-dimensional squeezing of the color-electric flux in the QCD vacuum. 
Hence, the QCD vacuum can be regarded as the dual version of the supercon- 
ductor. In this dual-superconductor picture for the QCD vacuum, the squeez- 
ing of the color-electric flux between quarks is realized by the dual Meissner 
effect as the result of condensation of color-magnetic monopoles. However, 
there are two large gaps between QCD and the dual superconductor picture. 

1. This picture is based on the abelian gauge theory, while QCD is a non- 
abelian gauge theory. 

2. The dual-superconductor scenario requires condensation of magnetic 
monopoles as the key concept, while QCD does not have such a monopole 
as the elementary degrees of freedom. 

As the connection between QCD and the dual superconductor scenario, 't 
Hooft proposed the concept of the abelian gauge fixingp, the partial gauge 
fixing which is defined so as to di- 
agonalize a suitable gauge-dependent 
variable The abelian gauge 

fixing reduces QCD into an abelian 
gauge theory, where the off-diagonal 
element of the gluon field behaves as 
a charged matter field (See Fig.l). 
As a remarkable fact in the abelian 
gauge, color-magnetic monopoles ap- 
pear as topological objects correspond- 
ing to the nontrivial homotopy group 
H 2 (SU(N c )/U(l) Nc - 1 ) = ZNc-i. Here, 
we assume abelian dominance, which 
means that the only abelian gauge 
fields with monopoles would be essen- 
tial for the description of the nonper- 
turbative QCD, and neglect the off- 
diagonal elements, which is called 
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Fig.l Dual Superconductor Picture 
in QCD 
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abelian projection. Thus, by the abelian gauge fixing and the abelian projec- 
tion, QCD is reduced into abelian projected (AP-)QCD, which is abelian gauge 
theory including monopoles. If the monopole condenses, the scenario of color 
confinement by the dual Meisser effect would be a realistic picture for confine- 
ment in QCD. In this paper, using the lattice QCD simulation, we study the 
key of the dual Higgs theory: abelian dominance (Sec. 2), the appearance of 
monopole (Sec. 3) and monopole condensation (Sec. 4). 



2 Abelian Dominance in the Maximally Abelian Gauge 

Abelian dominance on the confinement force has been investigated using the 
lattice QCD simulation in the maximally abelian (MA) gaugdra. In terms of 
the link variable U^(s) = U®(s) + ir a U^(s), the MA gauge fixing is defined 
by maximizing R = £ s>/ , tr(U^s)r 3 Ul(s)r 3 ) = E^{(^(*)) 2 + (^(*)) 2 " 
(^( s )) 2 — (U^(s)) 2 } through the gauge transformation. In this gauge, the 
off-diagonal components, and U^, are forced to be small, and therefore 
the QCD system seems describable by U(l)-like variables approximately. It 
is to be noted that the MA gauge is a sort of the abelian gauge, because 
^K s ) = Yla ± ^±m( s ) t 3^±m( s ) ^ s diagonalized there. (For the simple notation, 
we use U-^(s) = L/t(s— /i) in this paper). In this section, we study the origin of 
abelian dominance on the confinement force in the QCD vacuum, considering 
the relation to abelian dominance on the link variable. 

2.1 Abelian Dominance on Link Variable: Microscopic Abelian Dominance 

In the lattice formalism, the SU(2) link variable L^(s) is factorized as 

TT . . / cos0 M (s) -sM^sJe-^WVe^W \ _ 
U ^=[ S me^s)e^ easels) ){ e -, e>) ) = M ^W[s), 

where u M (s) and M^(s) correspond to the diagonal part and the off-diagonal 
part, respectively. For the residual U(l)3 gauge transformation, u^s) behaves 
as the abelian gauge field, while c fJi (s) behaves as the charged matter. Here, 
the range of angle variables are taken as < 6> M < ^, — ir < 8^, x M < ir. 

In order to investigate abelian dominance on the link variable in the MA 
gauge, we define "abelian projection rate" as RAbel(s,fj) = cos0 M (s)a. For 
instance, if cosO^ = 1, the SU(2) link variable is completely abelian as i/ M (s) = 

ift 3 I , while it is off-diagonal as Uuls) = =0,3,^, it 
e " J \e ^ ** J 

cos# M = O.We show in Fig. 2 the spatial distribution of the abelian projection 

rate i?Abei = cos# as an arrow (sin 6, cos 9). In the MA gauge, most of all 

SU(2) link variables become U(l)-like. We show also in Fig. 3 the probability 
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Fig. 2 Local abelian 
projection rate -f?Abd = 
cos at 0= 2.4 on 16 4 
lattice. Arrows denote 



(sin 0, cos 0) [0 < < §]. 



Let us consider the Wilson loop {W C [U ^{s)]) = (trll c U^s)) in the MA 
gauge. In calculating {Wc[U l _ l (s)]), the expectation value of e* x ^ s ) in M M (s) 

vanishes as (e'^^'jMA — f " dX/j,{s) exp{ix M (s)} = 0, when x^( s ) ' 1S assumed 
to be a random angle variable. Then, the off-diagonal factor U^(s) appearing in 
the Wilson loop WclU^s)] becomes a diagonal matrix, Z7^(s) = M /J (s)u Jlt (s) — > 
cosf? M (s)u AI (s), in the MA gauge. 

Then, for the J x J rectangular C, the Wilson loop Wc[U fl (s)} in the MA 
gauge is approximated as L L L 



(Wc[U„(s)]} = (trJJ ~ <JJcoBfl w («i) -tr 



MA 



(exp{> ln(cosf?^ i (si))})MA (Wc[%(s)])ma 




exp{L(ln(cos0 M (s)))MA} (W / c[«m( s )])ma, 



(1) 
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Fig. 3 The probability distribution 
f(-RAbei) of abelian projection rate 
RAbei{s,fi) at (3 = 2.4 on 16 4 lattice 
from 40 gauge configurations in the MA 
gauge (solid curve) and without gauge 
fixing (dashed curve). 




L = 2(I+J) 



Fig. 4 The comparison between the 
analytical estimation (straight line) 
and the lattice data (x) of W(f = 
(Wc[U])/(W c [u])ma with I, J > 2 as 
the function of the perimeter L = 2(7 + 
J) in the MA gauge at [3 = 2.4. 



where L = 2(1+ J) denotes the perimeter length and Wcfw^s)] = ^ T ^-i=i u in i s i) 
the abelian Wilson loop. Here, we have replaced X)i=i ln{cos(# Mi (s^)} by its av- 
erage L(ln{cos ^(s)})ma in a statistical sense. In this way, we derive a simple 
estimation as Wf = (Wc[U^{s)})/(W c [u^{s)})ma ^ exp{L(ln(cos^(s))) M A} 
for the contribution of the off-diagonal gluon element to the Wilson loop. From 
this analysis, the contribution of off-diagonal gluons to the Wilson loop is ex- 
pected to obey the perimeter law in the MA gauge for large loops, when the 
statistical treatment would be accurate. 

In the lattice QCD, we find that Wg s seems to obey the perimeter law for 
the Wilson loop with /, J > 2 in the MA gauge, as shown in Fig. 4. We find 
also that the lattice results on W^ s as the function of L is well reproduced 
by the above estimation with microscopic information on the diagonal factor 
cos ^(s) as (ln{cos M (s)}) M A ^ -0.082 for [3 = 2.4. 

Thus, the off-diagonal contribution Wg ff to the Wilson loop obeys the 
perimeter law in the MA gauge, and therefore the SU(2) string-tension <7su(2) 
becomes same as the abelian string-tension fTAbci, 

ct S u(2) = - lim ln(W Jx AU^s)}) 

I, J— >oo 1 J 

~ -2(ln{cos0u(s)}) M A lim - ^/ + er A bci /; ^>°° <7 A bei, (2) 

when the finite size effect on I and J is removed. Thus, abelian dominance for 
the string tension, crsu(2) — fAbci, can be proved in the MA gauge by replacing 
the off-diagonal angle variable Xm( s ) as a random variable. 
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3 Monopoles in QCD 



In this section, we derive the abelian gauge theory including monopoles in 
terms of the gauge connection, and study monopole properties! 

In the general system including singularities such as the Dirac string, the 
field strength is defined as G^ v = i([Z) M ,Z)„] — [d^,d v ]) with the covariant 
derivative = + ieA^. By the general gauge transformation with the 
gauge function fi, the field strength G^„ is transformed as 



^linear , ^bilinear i ^sing /o\ 

The last term G™ g remains only for the singular gauge transformation, and 

can provide the Dirac string. 

Here, in order to study appearance of the monopole, let us consider the 

/ e i<p cos- sin- \ 
SU(2) singular gauge-function — [ . fl 2 _..,„ 2 fl ) with polar angle 

\ — sm| e * v cos| J 

9 and azimuthal angle tp. Here, fl is multi-valued at the positive region of 
z-axis, 8 — 0, and the last term G™ s in Eq.(|^) provides the singular gauge 
configuration of the Dirac string in the abelian sector, as shown in Fig. 5. Ac- 
cordingly, the linear term Gjj£ ear and the second term G^ mcar in Eq. (||) include 
the monopole with the Dirac string and the anti-monopole, respectively. By the 
abelian projection — > = A^T 3 , the bilinear term, which is originated 
from the non-abelian nature, is dropped. The field-strength in the abelian 
projected QCD is derived as F M „ = d^A v — d v A^ + d u ]£ft, which breaks 

the abelian Bianchi identity. Thus, by the singular gauge transformation as f2, 
the monopole appears in the abelian sector of QCD. 

Next, we consider monopole properties in terms of the action. In the 
static frame of the monopole with the magnetic charge g, a spherical 'mag- 




Abelian Projected QCD 
Fig. 5 Appearance of monopoles in abelian projected QCD. 



netic field' is created around the monopole in the abelian sector of QCD 
as H(r) = with = eijkdjA 3 . Then, the monopole inevitably ac- 

companies a large fluctuation of the abelian gluon component A 3 ^ around 
it. For the abelian part S^bd = — | J d i x(d fJ ,Al - d v A^) 2 of the QCD ac- 
tion, the electro-magnetic energy created around the monopole is estimated 
as E(a) = f°° (Px^Hfr) 2 — ^— , where a is an ultraviolet cutoff like a lattice 

* / J a 2 \ J 87m 1 

mesh. As the "mesh" a goes to 0, the monopole provides a large fluctua- 
tion of SAbei, and hence the monopole seems difficult to appear if the abelian 
gauge theory is controlled by SAbei- This is the reason why QED does not 
have the point-like Dirac monopole. Then, why can the monopole appear in 
the abelian sector of QCD ? To answer it, let us consider the division of the 
total QCD action Sqcd into the abelian part S*Abci and the remaining part 
5 ff = Sqcd — SAbei, wn i cn is contribution from the off-diagonal gluon compo- 
nent. Unlike Sqcd and SAbei, '-'off is not positive definite and can take a neg- 
ative value in the Euclidean metric. Then, around the monopole, the abelian 
action SAbei should be partially canceled by the remaining off-diagonal contri- 
bution S ff such that the total QCD action Sqcd around the monopole does 
not become extremely large. Thus, we expect large off-diagonal gluon com- 
ponents around the monopole for cancellation with the large field fluctuation 
of the abelian part. Based on this analytical consideration, we study action 
densities around monopoles in the MA gauge using the lattice QCD. 

On the SU(2) lattice, we measure the action densities Ssu(2)> SAbei and 
S ff, which are the SU(2), the abelian and the off-diagonal parts, respectively. 
We show the probability distributions of Sgu(2), SAbei and S s in Fig. 6. (To 
be exact, S is the averaged value over the neighboring links around a dual 
linlM) For the total distribution on the whole lattice, most S g- are positive, 
and both SAbei and S Q ff tend to take smaller values than Sgu(2) owing to 




S(m) 

Fig. 6(a) Probability distributions P(S) 
of the Sgij(2) (dashed curve), SAbcl (solid 
curve) and S ff (dotted curve) in the MA 
gauge at /3 = 2.4 on 16 4 lattice. 
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Fig. 6(b) Probability distribution P k (S) 
around the monopole for S , gu(2)i ^Abel 
and S a ff in the MA gauge at /3 = 2.4 on 16 4 
lattice. 
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^su(2) = ^Abei+'Sofi- Around the monopole, however, the off-diagonal part S s 
of the action density tends to take a large negative value, and S a e strongly 
cancels with the large abelian action density S'Aboi so as to keep the total 
SU(2) action small. Thus, monopoles can appear in the abelian sector in QCD 
without large cost of the QCD action due to the strong cancellation between 
the abelian action SAbci and the off-diagonal part S Q s of the action. 

In conclusion, in the small distance, all SU(2) components contribute to 
the short-distance physics. In the dual-superconductor picture, the off-diagonal 
components are necessary to create the monopole in the abelian sector. In the 
large distance, however, the off-diagonal components are not necessary for the 
infrared physics, and only abelian gauge fields with monopoles describe the 
QCD vacuum. 



4 Monopole Condensation 



In the dual Higgs theory, the monopole is assumed 
to be condensed, which has been suggested by the 
formation of the global network of the monopole cur- 
rent in the lattice QCD. In this section, we consider 
monopole condensation in the QCD-vacuum at the in- 
frared scale. However, QCD is described by the gluon 
field not by the monopole current, and therefore it is 
difficult to clarify monopole condensation only with 
the lattice QCD simulation. To this end, we generate 
the monopole-current system on the lattice using a 
simple monopole-current action, and study monopole 
condensation and the role of the monopole to color 
confinement]. 

In general, monopole-current action S[fe M (s)] in- 
cludes the nonlocal Coulomb interaction as Sc = 
J d 4 xd 4 yk fJl (x)D(x — y)k ll (y) with the Coulomb prop- 
agator D{x). In the dual Higgs phase, however, 
the effective interaction between monopole currents 
would be short-range due t4hithe screening effect by 
the dual-EJiggs mechanism Ettl similar to the Debye 
screening!^. Then, the infrared partition functional 



of the monopole current k^x) 
written as 



— 2.3 



would be 




Fig. 7 Monopole current 
in the monopole-current 
system in R 3 at a fixed 
time. 



Z = y Dfc M cxp{— a J d 4 xtrk^ l (x)}S(d fl k ll ), 



(4) 
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where a is the energy per unit length of the monopole current. Here, a is an 
ultraviolet cutoff larger than the screening length. We put the system on the 
4-dimensional lattice with the lattice spacing a. Fig. 7 shows behavior of the 
monopole current in the monopole current system for a = 1.5, 1.8, 1.9. When 
energy a becomes smaller, the monopole current is generated more easily. As 
a = 1.5, the monopole current covers the whole lattice and creates the global 
network, which is similar to the monopole current in the QCD vacuum. 

Here, let us consider this behavior of monopole currents in the analytical 
way. In_the abelian gauge of QCD, the charge of the monopole is expected to 
be unitjfl. The partition function can be approximated as the single monopole- 
loop ensemble with the length L, Z — X)iP(^) e ~ aL , where L and p(L) are 
length of the monopole loop and its configuration number, respectively. The 
monopole current with the length L is regarded as the L step self-avoiding 
random walk, where 2d — 1 = 7 direction is possible in each step. Therefore, 
p(L) is roughly estimated as (2d — 1) L = 7 L , and the expectation value of the 
monopole-current length is found to be 



When energy a is larger than "entropy" ln(2c? — 1), the monopole-loop length 
is finite. However, when a is smaller than entropy, the monopole-loop length 
becomes infinite, which corresponds to monopole condensation in the current 



\n(2d — 1) ~ ln7 ~ 1.945, which corresponds to "entropy" of the self-avoiding 
random walk. Such a transition is quite similar to the Kosterlitz-Thouless 
transition in 2-dimensional superconductors, where vortex condensation plays 
an important role to the transition. 

Next, we study how these monopole currents contribute to color confine- 
ment properties! Quark confinement is characterized by the linear inter-quark 
potential, which can be obtained from area-law behavior of the Wilson loop, 
(W) — (Pexp(iej> A^dx^)). The expectation value of the Wilson loop (W) 
is shown in Fig. 8. The Wilson loop exhibits the area law and the linear con- 
finement potential: h^VF) decreases linearly with the quark loop area I x J. 
Quantitatively, the string tension is measured by the Creutz ratio, and we 
show in Fig. 9 x(3, 3) as a typical example. For the monopole condensed phase 
as a < a c , the string tension gets a finite value, while it vanishes for the 
non-condensed phase of monopoles as a > a c . Thus, the confinement phase 
directly corresponds to the monopole condensed phase, and therefore monopole 
condensation is considered as essence of the confinement mechanism. 




(5) 



representatio: 



Here, the critical value on monopole condensation is a c 
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Fig. 8 The Wilson loop (W{I X J)) in the 
monopole-current system. 
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Fig. 9 The Creutz ratio as the function of 
a in the multi-monopole system. 



We compare the lattice QCD with the monopole-current system in terms 
of monopole condensation and confinement properties. The lattice QCD simu- 
lation shows that the QCD vacuum in the MA gauge holds the global network 
of the monopole-current. Considering the similarity on the monopole clus- 
tering, the QCD vacuum can be regarded as the monopole condensed phase 
with a < a c in the monopole-current system, as shown in Fig. 7(c). Such 
an identification of the QCD vacuum with the monopole condensed phase is 
also suggested in terms of the confinement properties, because the confinement 
phase corresponds to the monopole condensed phase as shown in Fig. 93. Thus, 
monopole condensation is considered to occur at the infrared scale in the QCD 
vacuum and the confinement would be caused by monopole condensation there. 

One of authors (H.I.) is supported by Research Fellowships of the Japan 
Society for the Promotion of Science for Young Scientists. 
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